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All data are referenced to fuselage body axes according to right-handed sign 
conventions . 
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Chapter 1 


INTRODUCTION 


In September of 1981 the National Aeronautics and 
Space Administration at the Ames/Dryden Flight Research 
Facility began the investigation of the "Deep Stall" 
phenomenon using a modified SGS 1-36 sailplane. This 
investigation was directed at demonstrating the 
feasibility of unpowered Controlled Deep Stall flights at 
ang 1 es-of -at tack of between 30 and 75 degrees. The primary 
research objectives of the SGS 1-36 Controlled Deep Stall 
Sailplane Project was to control the aircraft in the Deep 
Stall regions by using large in-flight deflections of the 
a 1 1 mo veable horizontal stabilizer, and to assess the 
vehicle's longitudinal and lateral-directional dyn am i c 
stability and control characteristics. Initial qualitative 
flight tests were conducted on radio controlled mode I sail- 
planes having both low and high (T-tail) horizontal 
stabilizers. A quarter scale model of the modified SGS 1-36 
was later used in wind tunnel tests at NASA Langley 
Research Center to obtain the static stability and control 
aerodynamic characteristics of the aircraft. These data 
along with the estimated rotary derivatives were 
incorporated into a fixed based flight simulation to 
facilitate engineering evaluation, pilot familiarization, 
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and planning for the full scale flight tests. Due to 
limited flight operations support and available time, only 
a total of twenty manned flight tests were conducted using 
the modified SG5 1-36 at C.G. positions of 33.1% and 28.4% 
Mean Aerodynamic Cord (MAC). 

Since the early 1940's, several variations of the 
Deep Stall concept have been incorporated into "free 
flying" model airplanes to recover them rapidly and safely 
at the end of their flight by "popping" the horizontal tai 1 
using a timing device. This technique, called the 
"Detherma 1 i zer" (Reference 1) , first introduced in 1942 by 
Carl Goldberg, has since generated interest toward the 
possible applications in safe recovery of high aspect ratio 
r emo tely piloted vehicles fr om high altitudes through 
turbulent atmospheric layers. In addition, in recent years 
several modern transport aircraft have experienced 
accidents related to the Deep Stal 1 phenomenon either in 
service or during flight test programs. Despite the 
incorporation of computer-contro 1 1 ed stick pushers on many 
of these aircraft to prevent entering an ang 1 e-of -at tack 
region where a Deep Stall could occur, several "T" tail 
configured aircraft (such as the British Aerospace BAC 111, 
Boeing 727, Canadair Chal lenger CL-600, F 101, F 104, and 
several high performance sailplanes) have crashed when they 
entered a stable deep stall where the pilot was not able to 
recover by means of normal stall 'recovery techniques 
because of inadequate pitch control (Reference 2). 
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Given adequate pitch authority through the use of 
large in-flight horizontal tail deflection, an uncontrolled 
flight may then become cont ro 1 1 ab 1 e. This increase ~in pitch 
control authority can be used in assisting to recover from 
inadvertent deep stall, spins, or spiral dives in IFR 
conditions. Additionally, near vertical descents in areas 
of restricted lateral maneuvering and precision recovery of 
Remotely Piloted Veh i c 1 es (RPV ' s ) in conjunction with a 
retro-rocket landing system could be achieved using the 
controlled deep stall technique. 

This report describes the flight test procedure and 
discusses the preliminary analysis of the results obtained 
from twenty manned flights of the SGS 1-36 in the high 
angles-of-attack Deep Stal 1 region. A comparison of the 
flight determined stability and control derivatives and 
those of the wind tunnel and the estimated aerodynamic data 
is also presented. Furthermore, deep stall dynamic response 
of the SGS 1-36 is discussed briefly to explain s ome of the 
unexpected flight observations. 


Chapter 2 


VEHICLE DESCRIPTION AND MODIFICATIONS 

The NASA SGS 1-36 Control led Deep Stal 1 Vehicle is a 
mo dified version of the standard single place, T-tail 
design Schweizer 1-36 sailplane that is commercially used 
as an advanced trainer. It is made primarily of aluminum 
structure and skin except the rudder surface that is fabric 
covered. The modifications were made to the fixed horizon- 
tal stabilizer, control cables, canopy, cockpit, and the 
nose cone of the aircraft. 

The physical characteristics of the vehicle are 
presented in Tables 1 through 3. The aircraft moments of 
inertia listed were obtained experimentally by means of an 
"inertia swing". Measurements were made at 33.1% mean aero- 
d y n am ic cord and the resulting inertia were analytically 
corrected to 23.8% mean aerodynamic cord C.G. position for 
data reduction of flights at this forward center of gravi- 
ty. Figure lisa picture of the modified SGS 1-36. Figure 
2 is a three-view drawing of the modified vehicle. 

Aerodynamic control is accomplished by means of 
conventional manual control system. The pilot is provided a 
center stick for elevator and aileron control, pedals for 
rudder control, and a side control lever to deflect the 


horizontal stabilizer. In addition, a spoiler control 
handle is provided for speed control and rapid descent. 
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The fixed, single control, standard T-tail was modi- 
fied into an all moveable dual control horizontal tail 
configuration. The new configuration allows the horizontal 
stabilizer to be deflected fromO to 73 degrees trailing 
edge up (Figure 3) through the use of the horizontal stabi- 
lizer control lever mounted on the port side of the cockpit 
(Figure 4). Only partial elevator travel authority is 
retained throughout the range of stabilizer deflection. As 
the stabilizer setting is increased to 73 degrees trailing 
edge up, the elevator travel authority is reduced from ful 1 
-24 to +17 degrees deflection to only -11 to -24 degrees. 
Furthermore, two six-pound counter weights were added to 
the stabilizer to reduce the risk of in-flight flutter of 
the tail section (Figure 3). 

During the original modifications to the horizontal 
stabilizer, cables and pulleys were employed to connect the 
stabilizer control lever to the stabilizer, but due to 
excessive cable elongation and friction under load, the 
cable and pulley system was replaced by a bellcrank and 
pushrod system. 

In order to assure safe and quick pilot exit from 
the cockpit in case of an uncontrolled flight, several 
modifications were made to the cockpit and canopy. An 
"Egress" pneumatic escape system was added to the the 
cockpit that would enable the pilot to eject the canopy and 


release his shoulder harness and seatbelt simultaneously, 
by pulling a single escape handle located to his right on 
the instrumentation panel (Figure 5). 

The original nose cone was modified to acconrmodate a 
special made nose boom. Figure 6 contains a picture of the 
nose cone with the boom installed. 


Chapter 3 


INSTRUMENTATION 


In flight data were obtained using the onboard data 
acquisition system and were transmitted and recorded digi- 
tally on magnetic tapes at the ground station. This system, 
a 28-channel 10-bit pulse code modulation tel erne ter y 
system, consisted of a 3-axis rate gyro, a vertical gyro, 
3-axis linear accelerometers, and control position 
transmitters (CPT's) to measure angular rates, Euler 
angles, linear accelerations, and control surface positions 
respectively. A1 1 of the system components are mounted on a 
single a 1 uminum platform for easy installation and remova 1 . 
Figure 7 depicts the data acquisition platform onboard the 
a i rcraf t . 

Ang 1 e-of -at tack, ang 1 e - o f - s i de s 1 i p , dynamic pres- 
sure, and static pressure were measured using a special 
made nose boom. This nose boom contained a self-al ignlng 
pitot-static tube, ang 1 e-o f - a 1 1 ack and ang 1 e - o f - s i de s 1 i p 
vanes and was wired to the data acquisition package onboard 
the aircraft (Figure 8). Corrections were made to the 
a n g 1 e-of-attack and angle-of-sideslip data for boom 
position and alignment. The control surface position, 
angular rate, and translational acceleration data were 
taken at 110, 220, and 440 samples per second (SPS) 
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respectively. For post-flight data analysis the sample rate 
used was 50 samples per second for longitudinal, and 25 SPS 
for lateral-directional data. In addition to the basic 
flight instrumentation six flutter accelerometers were 
installed externally on the vertical and horizontal tail 
surfaces for initial ground and flight flutter clearance 
tests. 

Table 4 includes the instrumentation parameters and 
their corresponding range, resolution, and accuracy. 


Chapter 4 


FLIGHT TEST PROCEDURES AND MANEUVERS 


A total of twenty manned flight tests were conducted 
at C.G. positions of 33.1% and 28.4% Mean Aerodynamic Cord 
using the modified NASA SGS 1-36 vehicle. Prior to the 
first deep stall flight the C.G. position of 33% MAC wa s 
selected to assure that the aircraft will have adequate 
longitudinal stability in the conventional and deep stall 
flight regime (The standard SGS 1-36 is certified by the 
Federal Aviation Atfcnin i s t rat i on for flight at C.G. range of 
20% to 40% MAC). It was also decided to fly the aircraft at 
the C.G. position of 28% MAC to increase the longitudinal 
stability and obtain trim data at this center of gravity 
configuration. Weight and balance measurements were made 
to verify the location of the vehicle's longitudinal center 
of gravity prior to initial flights at the two desired C.G. 
locations. Moments of inertia were determined 
experimental ly at 33.1% mean aerodynamic cord by means of 
an "inertia swing" and corrected analytically for 28.4% 
mean aerodynamic cord center of gravity position. 

Figure 11 is the complete flight schedule matrix 
i nd icat i ng the number of flights and their respective C.G., 
stabilator setting, and control stick position. 
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As a glider, the SGS 1-36 was towed to an approx i - 
ma te altitude of 3353 meters (11,000 feet) me an sea level 
(MSL) and then released for an unpowered flig’ht. Upon 
release, the pilot would maneuver the aircraft into a 
steady-state flight condition. The airspeed was then 
decreased to the stall speed of about 36 Knots using the 
center control stick. Just prior to complete stall, the 
stabilizer control lever was moved back quickly but gently 
to a preset position of 40, 50, or 60 degrees stabilizer 
setting depending on the required ang 1 e-of -at tack. This 
action caused the angle-of-attack of the aircraft to 
increase rapidly through the "transition region" of between 
15 to 30 degrees to the desired deep stall angle-of-attack. 
This so-called transition region, is characterized by large 
asymmetry in the lateral-directional forces and moments 
(more pronounced for high aspect ratio aircraft) that would 
cause the SGS 1-36 to enter an undesirable spin or spiral 
dive. 

Once the aircraft was in the deep stall region, the 
pilot stabilized the aircraft on a trim point for about 
five seconds and then performed elevator, rudder, and aile- 
ron doublet maneuvers. As the next step, handling quality 
evaluation maneuvers were performed using the elevator, 
rudder, and ailerons. At an approxi ma te altitude of 2393 
meters (7850 feet) MSL or 1676 meters (5500 feet) above 
ground level (AGL) recovery was initiated by using the 
stabilator to lower the angle-of-attack to the transition 



area. Although the actual ang 1 e-of -at tack was not displayed 
in the cockpit, the transition region was sensed by the 
pilot from the mild tail buffet due to the wing Wake. Upon 
encountering the tail buffet, the stabilator was rapidly 
moved to zero stabilator deflection position. During reco- 
very 91 to 152 meters (300 to 500 feet) of altitude were 
lost while the airspeed typically increased from 38 to 55 
Knots . 

The concept of the Controlled Deep Stall flight is 
illustrated in Figure 9. Figure 10 depicts the SG5 1-36 in 
deep stall flight at 50 degrees ang 1 e-of -at tack. 

During initial flight test planning, provisions were 
made to fly the aircraft at 30, 40, 50, 60, and 70 degrees 
trailing edge up stabilizer setting with the C.G. at two 
different positions. Initial ly, the wind tunnel data had 
indicated that with the C.G. at about 33% MAC there would 
be a one-to-one correspondence between the stabilizer 
setting and the ang 1 e-of -at tack. However, during the first 
and second deep stall flights with the horizontal stabili- 
zer at 60 degrees and the center of gravity at 33.1% MAC, 
the average ang 1 e-of -at tack was recorded to be 72 degrees. 
In addition, some unexpected 1 ater a 1 -d i rect i ona 1 , apparent 
Dutch-Roll, oscillations were observed. Based on this 
observation the remaining flights with higher stabilizer 
setting (70 degrees, with forward and aft center stick 
positions) were deemed too risky and were consequently 
cancelled. With the stabilizer setting at 40 degrees 
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trailing edge up, C.G. position of 28.3% MAC, and elevator 
control stick fully aft the pilot was not able to maintain 
the aircraft in the deep stall region for any s'ustained 
period of ti me . In this flight regime the aircraft had a 
tendency to enter and r ema in in the transition region. 
Hence, it was decided not to perform the remaining one 
flight in the s ame configuration with the stick fully 
f o rwa r d . 

During the first deep stall flight of the SGS 1-36 
with the a ng 1 e - o f - a t t a c k of about 72 degrees, two 
unexpected phenomena were encountered. Firstly, as soon as 
the aircraft entered the deep stall flight, the pilot 
experienced an unusually large, unstable lateral control 
stick hinge moment. He estimated the force on the center 
stick to be approximately 25 pounds and indicated that he 
had difficulty keeping the stick centered using his right 
hand (he used his left hand to keep the stabilator control 
lever in the deflected position) since the stick had a 
tendency to deflect either to the right or to the left. 
Secondly, as men t ioned earlier, in the deep stall flight 
the aircraft exhibited an unstable 1 atera 1 -d i rect i ona 1 
oscillation. During the flights that followed it was 
discovered that the magnitude of the unstable control stick 
hinge moments and the lateral -directional instabilities 
were functions of the ang 1 e-of -at tack. As the angle-of- 
attack increased so did the magnitudes of the unstable 
stick hinge moment and the apparent divergent Dutch-Rol 1 
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oscillation to a constant limit cycle. This phenomena will 
be discussed in more detail in the following sections. 

Due to rapid aircraft descent in the de.ep stall 
region, up to 1372 meters per minute (4500 feet per 
minute), the pilot had only about 60 to 100 seconds to 
complete his maneuvers. Because of the limited time availa- 
ble for obtaining flight data and rapidly changing flight 
conditions, there was only one opportunity to perform each 
maneuver. Thus maneuvers were practiced on the flight simu- 
lator prior to each flight. Post flight analysis of these 
maneuvers indicated that a doublet or pulse, fol lowed by 
two to five seconds in which the pilot made no input, was 
most effective in providing satisfactory stability data. It 
was also discovered that the range of the ang 1 e-of -s ides 1 i p 
calibration was not sufficiently large, which prevented 
the magnitude of raw s i de - s 1 i p-ang 1 e s greater than 35 
degrees from being recorded. 
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WIND TUNNEL TESTS 

A low speed wind tunnel investigation was conducted 
at NASA Langley Research Center's 30 x 60 wind tunnel to 
determine the static aerodynamic characteristics of a 
quarter scale model of the modified SGS 1-36 sailplane 
(Reference 3). For economical reasons, the model used was a 
low-cost "free-f light" model which limited the wind 
tunnel's dynamic pressure and hence the Reynolds number. 
Longitudinal and 1 atera 1 -d i rect i ona 1 force tests were 
conducted over an ang 1 e-of -at tack range of 0 to 90 degrees 
and ang 1 e-of - s i des 1 i p range of +/-10 degrees. Control 
effectiveness was determined for horizontal stabilizer 
deflections of 0 to 70 degrees trailing edge up, elevator 
deflections from +20 to -25 degrees, aileron deflections of 
32 degrees up and 12 degrees down, and full rudder deflec- 
tion of +/- 27 degrees. The wind tunnel speed was measured 
at 12.19 meters per second (40 feet per second) which 
corresponds to a dynamic pressure of about 85.9 Newton per 
square meter (1.8 pounds per square foot) and Reynolds 
n umb er of about 1.8 x 10^, based on a chord length of .24 

meters (0.8 feet). This low Reynolds number is general ly 
considered to be on the high side of the sub-critical 
range. Sub-critical Reynolds number data is known to 
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greatly underpredict the effectiveness of any trailing edge 
control surfaces. 

During the initial wind tunnel test a significant 
change in the wing dihedral was observed due to extreme 
wing flexibility. Therefore, a limited second series of 
tests were conducted with the wings restrained from bending 
so that a comparison could be made between the two series 
of tests. The final test results were a combination of the 
two series of tests. All the wind tunnel data were provided 
in the aircraft stability axis and were later transformed 
into body axis for the final comparison with the flight 
test results. 

With regards to the validity of using low Reynolds 
number wind tunnel data for application to the full scale 
sailplane, Figure 12 of Re ference k indicates that the 
effects of Reynolds nianber on the longitudinal characteris- 
tics of the aircraft are minimal above 30 degrees angle-of 
-attack. Sincethe the aerodynamic characteristics of the 
SGS 1-36 in the range of 30 to 75 degrees ang 1 e - o f - a 1 1 ac k 
were of the prime interest, it was assumed that the final 
wind tunnel results were reasonable representation of the 
conditions that would exist at the full scale Reynolds 
number range of about 1.0 x lO^to 3.0 xlO^. The values of 
obtained from wind tunnel tests were corrected for high 
Reynolds nionber of the full scale aircraft based on results 
of Figure 12. 
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COMPUTED AERODYNAMIC STABILITY DERIVATIVES 

In the conduct of the SGS 1-36 Controlled Deep Stall 
Project, a complete set of st-atic and dynami c( r otar y ) 
aerodynamic stability derivatives were required for the 
ang 1 e - o f - a 1 1 ack range of -5 to 90 degrees, so that a 
comprehensive comparison could be made between flight test 
and predicted data. The wind tunnel testing of the SGS 1-36 
quarter scale mode 1 fulfilled part of the requir emen t by 
providing a set of static aerodynamic derivatives. In order 
to obtain a complete set of rotary derivatives, a decision 
had to be made as to which technique should be employed to 
predict these derivatives at such wide range of angles-of- 
attack. Theory had been we 1 1 established to estimate the 
rotary derivatives at low ang 1 es -of -at t ack (-5 to 15 
degrees), but there were no techniques available for deter- 
mining these derivatives at high ang 1 es-of-at tack range of 
20 to 90 degrees. Based on this fact, a decision was made 
to use a combination of computer programs and computational 
methodsto obtain the rotary derivatives in the low angles- 
of-attack, and to develop some new techniques to compute 
these derivatives in the high ang 1 es-of-attack region. 

Vortex-Lattice program of Reference 5 was used to 
compute Clp, Cn r , Cm^, CNq, and CY r in the low angle-of- 
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attack region of -5 to +15 degrees. This program, an 
extension of the finite step lifting-line method, assumes 
steady state, irrotational, inviscid, incompressible, 
attached flow and uses Prandt 1 -GI aver t similarity rule to 
model the compressibility effects. In the past, the Vortex- 
Lattice program has been used extensively at NASA 
Arne s/Dryden and the results indicated a good correlation 
between derivatives obtained using this program and that of 
the flight tests. Therefore, it was assumed that acceptable 
results could be obtained using the Vortex-Lattice program. 

The computational method used to determine the 
values of Cnp and Cl r at angles-of- attack of between - 5 and 
+15 degrees were based on the well established, empirical 
procedures of reference 6. 

To obtain results in the high angle-of-attack range 
of 20 to 90 degrees, an equation was developed to compute 
the values of Clp assuming drag as the only source of 
aerodynamic force and moment (Appendix C). No other 
analytical equations were possible to develop without some 
gross oversimplifications that would have rendered them 
useless. Therefore, esti ma tions of the r ema ining rotary 
derivatives were made purely based on trends exhibited by 
the available static and dynamic derivatives in this high 
angle-of-attack region. 

In the deep stall flight regime, the parameter CN^ 
was predicted to behave similar to the calculated Clp which 
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is also a function of the wing's C^. Given that Cm^ is a 
function of moment arm l t (the distance between the wing 
and the tail quarter cord ), and that this moment arm 
decreases as a function of COS(a) as the ang 1 e-of -at tack 
increases, it was assumed that the magnitude of Cm^ would 
also decrease as a function of COS(a) in the range of 20 to 
90 degrees. 

The value of CY r was selected to be constant through 
-5 to 90 degrees ang 1 e-of -at tack since the Vortex-Lattice 
computation of this parameter at low angles-of-attack 
region had resulted in a constant value. 

The parameters Cnp and Cl r were presumed to decrease 
with ang 1 e-of -at tack similar to the decrease in Clg a , 
aileron control effectiveness (a wing parameter), once past 
the stall ang 1 e-of -at tack. Similarly, Cn f was assumed to 
decrease with ang 1 e-of -at tack just as Cng r , rudder control 
effectiveness (vertical tail parameter), decreased above 20 
degrees ang 1 e-of -at tack. 

It is important to notice that the damping 
derivatives obtained using the Vortex- 1 att ice program were 
pure rotary derivatives without the inclusion of the 
t rans 1 at i ona 1 acceleration derivatives. In general, the 
d amp ing derivatives obtained fr om flight or oscillatory 
wind tunnel tests are a combination of the rotary and 
translational acceleration derivatives. As an example the 
two parameters Cm^ and Cm£ could not be individually 
determined from flight data, therefore the combination of 
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the two is commonly referred to as Cm^ , its rotary 
derivative name. In this report all flight derivatives are 
called only by their rotary derivative name. 

Ordinarily, during the damping derivatives calcula- 
tions, the translational acceleration derivatives are 
neglected because of their relatively small magnitude and 
lack of readily available techniques to compute them. 


However , s i nee translational acceleration derivatives Qr> £ 

and Cn^ often have significant magnitude, it is possible 

for the predictions of Cm_+Cm» and Cn_+Cni to be 1 ower than 

M a r p 

that of the actual flight due to under esti ma t i o n of Cm £ 


and Cn* . In predicting the values of Cm_+Cm.and Cn_+Cn«, 

p s a r p 

the Vortex-Lattice program was used to determine the 


magnitude of On^ and Cn r . Based on established guide lines, 
the magn i tudes of Cm^ and Cn^ were set equal to 4 0% of the 
values obtained for Qn^ and Cn f respectively. 

A complete set of wind tunnel determined static 


forces and moments as well as the predicted rotary aerody- 
namic stability derivatives versus ang 1 e-of -at tack are 
presented in the appendix A. 
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DATA ANALYSIS 

To determine the longitudinal and lateral- 
directional aerodynamic stability derivatives from the 
recorded flight test time histories a digital computer 
program was used. This Fortran program, the Modified 
Max imum Likelihood Estimation (MMLE) of Reference 7 is a 
generalized dynamic analysis program which includes 
specific provisions applicable to aircraft dynamic 
analysis. The MMLE program employs a modified Newton- 
Raphson iteration technique to minimize the error between 
the various flight recorded time responses and the 
corresponding responses of the mathematical model of the 
aircraft (Appendix B). As an example, the flight recorded 
time histories of Flight 7 are presented in Figure 13. 

The aircraft mathematical model is represented by 
two independent sets of three degree-of -f reedom linearized 
differential equations of motion with the stability deriva- 
tives as the unknown coefficients. These two sets of 
equations are known as the longitudinal and lateral- 
directional equations of motion. The longitudinal equations 
are excited by flight measured elevator deflection and the 
1 a te r a 1 - d i rec t i ona 1 equations by flight measured aileron 
and rudder deflections. The MMLE program uses the predicted 
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static and rotary derivative set as the starting 
coefficient values to solve either the longitudinal or 
1 ater a 1 -d i rect i ona 1 equations of motion and produce a 
system response time history. The resulting time responses 
and the corresponding flight time histories are then 
compared to determine the difference between them that 
constitutes the error. This error is represented by the 
error vector, 

AY= [AP*Ar»AP»Ar>A/3> Aa y 3 T 

where the objective is then to minimize AY in an optimal 

fashion using the cost functional 

t 

3 = /( AY)' W AY dt 

where W is a weight matrix reflecting the relative 
confidence in the inst rument measurement s. To accomp 1 i sh 
this task the Newton-Raphson iteration method was selected 
and modified to provide successful minimization. Once the 
error has been minimized the corresponding values of the 
coefficients are selected as the best estimation of the 
aircraft's aerodynamic stability and control derivatives. 
Figures 14 and 15 represent examples of MMLE produced and 
flight determined 1 ater a 1 -d i rect i ona 1 and longitudinal time 
histories for Flights 7 and 13. 

A measure of the accuracy of these estimated deriva- 
tives is provided in the form of the Cramfer-Rao bound. The 
smaller the Cramfer-Rao bound, the more confidence that can 
be placed in the estimated values of the derivatives. An 
evaluation of the use and accuracy of the Cramfer-Rao bound 
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is given in Reference 8. Ordinarily the Cramfcr-Rao bound 
produced by MMLE tends to indicate a high degree of 
accuracy for the estimated derivatives, where in' reality 
the accuracies are lower due to aerodynamic modeling errors 
(i.e. no consideration for non-linearities), poorly 
conditioned maneuvers, presence of state noise, and 
inevitable sensor errors (despite provisions within the 
MMLE to account for some of these modeling errors). To 
compensate for these anomalies the resulting Cramtr-Rao 
bounds were mu ltiplied by a factor of 3 based on previous 
experience that have shown this to produce a more realistic 
representation of the error bound . 

The MMLE program contains several options that 
facilitates the convergence of the initial values to a 
final derivative set. For stable systems, the program 
normal ly converges in 6 or 7 iterations, but in some 
instances where the dynamic response may be unstable or. 
non-linear, the program would not easily converge without 
the use of the "A priori" option. This option, ai lowed the 
starting set of derivatives to be weighted, which tended to 
hold the derivatives near their starting values if no 
information about them was contained in the maneuver. Early 
in the flight program, wind tunnel and the computed 
predictions were used as the starting values. However, as 
different trends in the data developed, previously obtained 
f 1 i ght -determi ned derivatives were used. Since a high "a 
priori" Value during complete iteration process tended to 
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produce unreal istic results, the weighing was removed after 
4 iterations so that the derivatives could change freely in 
the remaining 4 or 5 iterations as they converged to a 
final set of values. Because in the high ang 1 es-of-attack 
regions a small banking of the aircraft translated into 
some un r ea 1 i s t i ca 1 y high sideslip angles, the flight 
recorded ang 1 e-of -s ides 1 ip was multiplied by the cosine of 
ang 1 e-of -at tack to arrive at the corrected sideslip angle. 

To gain insight into the 1 a t e r a 1 - d i r e c t i on a 1 
instabilities observed during flights at high angle s-of- 
attack, the DIGIKON computer program of Reference 9 was 
used to conduct a simple open loop analysis of the aircraft 
dynamics using the results from Flight 7. DIGIKON 
represents an advanced state-of-the-art tool for modeling 
and analyzing digital and continous flight control systems. 
It uses the modern state-space approach for modeling the 
aircraft and its various flight control systems. To 
facilitate the aircraft modeling in the state-space form, 
an interactive interface program between the user and 
DIGIKON was created. This interface program, ACST, accepts 
the non-dimensional aerodynamic stability and control 
derivatives produced by MMLE, along with their correspon- 
ding flight condition and mass properties, to generate the 
d imens iona 1 i zed, linear, differential equations of motion 
of the aircraft in the state-space format, as expected by 
DIGIKON. Once the aircraft is modeled, DIGIKON al lows the 
user to interactively perform time domain or frequency 
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doma in analysis for a given input (control surface 
deflection). For this report only two S-plane Root-Locus 
plots were generated. One plot represents the aircraft 
response if the pilot (assumed to behave as a pure gain) 
tries to control the bank angle using his ailerons (Figure 
19), and the second plot, if he controls the sideslip angle 
using his rudder (Figure 20). 
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RESULTS AND DISCUSSION 

In the past, dynamic response and handling qualities 
of aircraft at very high ang 1 es-of -at tack have seldom been 
investigated. The NASA SGS 1-36 Controlled Deep Stall pro- 
ject has, for the first time, provided limited data in this 
extreme ang 1 es -of -at tack region for a full-scale, piloted 
aircraft. The analysis presented in this section is the 
first qualitative investigation of the phenomenon at high 
Reynolds numbers. Hence, the preliminary nature of the 
presented analysis is herein emphas i zed. 

Flight 7 has been selected as a case example to 
represent the results from the ma jority of the 20 flights 
that were conducted in the deep stal 1 region. For this 
reason, the data obtained fr om flight 7 is analyzed in 
detail in order to provide s ome explanation for the 
observed flight characteristics of the aircraft in this 
f 1 ight regime. 

Figure 13 is the presentation of the complete deep 
stall t ime histories for flight 7. An inspection of these 
time responses reveals the ease with which the pilot was 
able to enter and exit the deep stall flight. As the air- 
craft enters into the deep stall (the average angle— of 
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attack of about 65 degrees), the dynamic pressure decreases 
to about 239 N/m^ (5 psf). This low dynamic pressure 
eliminates any risk of possible structural damage to the 
airframe. The vertical descent rate increases to a value of 
1219 meters per minute or 4500 Ft/Min ( plotted data are 
truncated above 1219 m/min or 4000 Ft/min of vertical 
velocity due to sensor limitations). The elevator 
effectiveness is shown to diminish drastically at this high 
angle- of-attack, but the lack of any significant longitu- 
dinal oscillation indicates a damped pitch characteristic. 
One of the most important observations to be made is the 
1 ater a 1 -d i rec t i ona 1 oscillation that is signified by the 
bank and sideslip angles oscillations at the rate of 1.39 
rad/sec. Initially, the magn itude of the oscillations 
increases to a constant limit cycle as the aircraft enters 
the deep stall, with the ang 1 e - o f - a 1 1 ack of about 65 
degree, but decreases as the pilot decreases the angle-of- 
attack to about 50 degrees using his stabilator. The 
oscillations increase again as the ang 1 e-of -at tack is 
increased back to 65 degrees. This phenomenon indicates a 
direct correlation between the ang I e-of -at tack and the 
magnitude of the lateral-directional oscillation for 
flights in this high ang 1 e-of -at tack regime. 

The attitude ti me history plot points out a mo s t 1 y 
level aircraft all through the deep stall portion of the 
flight. The pitch rate time history denotes a small, higher 
frequency oscillation ( 3.14 rad/sec), but lower magnitude 
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(average of 14 deg/sec) than that of the rol 1 and yaw rates 
(1.4 rad/sec with average magnitude of 45 and 17 deg/sec 
respectively). The maximum normal acceleration of "about 1.9 
G is encountered only during the recovery phase, where as 
during the deep stall flight the normal acceleration 
r ema i n s a t 1 G. 

The apparent independence of the bank and sideslip 
angles from the aileron and rudder deflections connotes a 
very low aileron and rudder effectiveness at extreme angles 
-of-attack. As indicated in the previous sections the 
plotted sideslip angle ti me history is truncated due to 
limited calibration range at high angles-of-attack. 

Figure 14 is the MMLE computed 1 atera 1 -d i rect i ona 1 
time histories for a given maneuver of Flight 7 as compared 
to the flight determined time histories of the same flight. 
Three points could be made regarding these plots. First, 
the aileron and rudder doublet maneuvers are quite good 
despite the uncontrolled 1 atera 1 -d i rect i ona 1 oscillations. 
Second, the MMLE computed time responses match that of the 
flight quite we 1 1 . But the ma tches during the aileron and 
rudder deflection time intervals are not as good. This 
indicates that the final MMLE output may not provide a 
realistic control effectiveness values for this maneuver. 
Third, since the sideslip angle was truncated due to 
calibration limitations, a low weight was imposed on this 
parameter during the MMLE iterations to reduce its ill 
effects on the estimated derivatives. It should be noted 
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that the lateral acceleration data provide the same 
information as that of sideslip angle. Therefore, MMLE is 
not deprived of any significant information. The ’lateral- 
directional time hi story matches presented in Figure 14 are 
typical of the other matches produced by MMLE for the 
remaining 19 flights at high angl es-of-attack. 

The longitudinal flight determined and MMLE computed 
time histories of Flight 13 are compared in Figure 15. This 
figure is an example of a typical MMLE match for 
conventional low ang le-of -at tack flight (about 2.5 
degrees). The ang 1 e-of-attack time response, following the 
elevator doublet maneuvers, implies an over damped 
longitudinal dynamic characteristic for the aircraft in 
this flight regime. 

In Figure 16, the flight envelope of the SGS 1-36 is 
presented in the form of flight determined trim data 
(elevator deflection versus ang 1 e-of -a t tack) , at C.G. 
positions of 33.1% and 28.4% MAC, and four stabilator 
settings of 0, 40, 50, and 60 degrees. In addition, the 

normal and controlled deep stall flight regime as well as 
transition and uncontrolled deep stall flight regions are 
clearly marked by dashed line. However, it should be noted 
that the dashed lines only approximate the border between 
the different regions. 

One of the significant characteristics of the trim 
plots is that in the normal flight regime the slope of the 
elevator deflection versus ang I e-of -at tack decreases as the 
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center of gravity is shifted from 28.4% to 33.1% MAC, in a 
conventional manner. Where as, in the deep stall region no 
significant change of slope is noticed. Instead, the 
plotted trim data exhibit a shift to the right foil owing 
the C.G. shift to the aft. 

The final MMLE estimated and predicted longitudinal 
stability and control derivatives of the NASA SGS 1-36 
sailplane are presented in Figure 17. The lateral- 
directional stability derivatives are plotted in Figure 18. 
The presented MMLE estimates, consist of results obtained 
from the 20 flights of the vehicle at the two C.G. 
positions of 28.4% and 33.1% MAC, which have been corrected 
to 40% MAC for comparison to the wind tunnel and the 
estimated aerodynamic derivatives. Those MMLE estimates 
with high Cramfer-Rao bound are el iminated from the plots 
according to the tolerances set forth in Table 5. 

The flight determined values of the Cm a and Cm^ 
(Figure 17) disclose a more stable longitudinal 
characteristic for the aircraft than that of the 
predicted, in the deep stall flight regime. On the other 
hand, the conventional flight regime comparison indicates 
that the damping derivative Cm^ , was not predicted with 
accuracy using the Vortex-Lattice program. The plot of CN 
versus the angle of attack reveals a higher than predicted 
normal force in the deep stal 1 region. As denoted by Figure 
17, the elevator effectiveness decreases sharply as the 
ang 1 e-of -attack is increased to 72 degrees, but the 
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negative values of CNg e might be due to the fact that 
elevator effect could not be modeled separately, and its 
value is a combination of elevator and stibi 1 ator 
effectiveness (stabilator would deflect about 2 degrees due 
to aerodynamic moment produced by elevator deflection). 

The 1 ater a 1 -d i rect i ona 1 static and dynamic stability 
and control characteristics of the NASA SGS 1-36 sailplane 
as estimated by MMLE are compared to that of the predic- 
tions in Figure 18. The parameter CY^ , displays lower than 
predicted magnitudes in the deep stall as well as conven- 
tional flight regi me . In the deep stall this par ame ter 
displaysa magnitude approximately equal to that of the 
normal flight regime. Since in the conventional flight 
regime the bulk of the side force due to side-sl ip-angle is 
created by the vertical stabilizer, and that in the deep 


stall the air flow over the vertical stabilizer is reduced 


markedly, it would therefore be reasonable to assume that 
the bulk of the side force in the deep stall flight regime 
is created by the fuselage. 

It is important to notice that the parameters CYg r 
and Clg r do not contain as many flight determined data 
points as the rest of the parameters. This occurrence can 
be attributed to the lack of information contained in the 
flight data about these parameters which when supplied to 
MMLE results in high Cramfer-Rao bound. A Cramfer-Rao bound 
greater than the predetermined tolerances (Table 5), will 



31 


result in the omission of s ome of those derivatives fr om 
the presented plots. 

Although the quarter scaled wind tunnel model used 
the same amount of differential aileron deflection as the 
full scaled aircraft, the examination of CYg a at low 
ang 1 es— of— at tack points out to higher wind tunnel predica- 
tion than the MMLE produced estimates (differential aileron 
deflection can reduce or eliminate the ma gnitude of this 
parameter). The parameter CYg a is, also, a weak parameter, 
that is to say, the deep stall and 1 ow angles values for 
this par ame ter could not accurately be esti ma ted using 
MMLE. 


The undesirable, lower than predicted values for 
Cl p in the deep stall region, is perhaps due to the 
Reynolds number difference between the wind tunnel model 
and the actual aircraft. A combination of low MMLE esti- 
mated values of this effective dihedral parameter and Cn^ , 
a primary directional stability parameter, may be one the 
major contributors to the 1 atera 1 -d i r ect i ona 1 instability 
observed in the high ang 1 es -of -at tack flights. The low 
angles estimation of Cn^ corresponds well with the predic- 
tion, except the negative value of one of the ma n e u v e r s , 
that should be regarded as unrealistic. 

The plot of Cl^versus ang 1 e-of -at tack , follows the 
wind tunnel predicted trend quit accurately, except with 
lower magnitudes in the conventional flight regime. The 
value of this parameter approaches zero as the angle-of- 
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attack is increased to 72 degrees, which indicates lack of 
any significant aileron effectiveness in the extreme 
angles-of-attack, as expected. 

The tendency of the Cng a to change from a proverse, 
favorable yaw characteristic at low angles-of-attack to 
adverse undesirable yaw character i s t i cs in the high angles- 
of-attack region is emphasized in Figure 18. The favorable, 
proverse, yaw characteristics of the aircraft in the 
convent i ona If 1 i ghtreg ime is, primarily, due to the large 
differential aileron deflection of this high aspect ratio 
vehicle. This desirable handling quality diminishes rapidly 
at higher angles when the completely separated flow 
eliminates any aileron effectiveness. 

A higher than expected rudder effectiveness is 
displayed at lower angles, and with lesser degree in the 
higher angles-of-attack by the plot of Cng r . The flight 
determined values of Cl p in the conventional flight regime 
nicely fol lows the trend set forth by the prediction but 
with lower magnitude. On the other hand, as the aircraft 
enters the deep stall region the value of this primary roll 
damping derivative rapidly approaches zero. The same 
observation holds true for Cl r , the primary yaw damping 
derivative. A combination of near zero value of Cl p and Cl r 
is the principal cause of the 1 atera 1 -d i rect i ona 1 oscilla- 
tion at angles-of-attack higher than 60 degrees. Although, 
the low angles MMLE estimations of Cn p exhibits a definite, 
precise trend, the larger than predicted negative values 
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for this parameter at low angles-of-attack, is an 
undesirable handling quality for the aircraft. 

To gain further insight into the 1 atera 1 -d i rect i ona 1 

dynamic response of the 5GS 1-36 in the deep stall region, 

« transfer 

5r 

functions for Flight 7 are presented in Figures 19 and 20. 
Because most pilots use visual cues to fly aircraft, and 
that the bank and sideslip angles are two of the strongest 
visual cues, $/g a and /3/^ transfer functions were selected 
for the open loop analysis of the vehicle, to obtain the 
closed loop dynamic characteristics of the aircraft (the 
pilot is assumed to be pure gain). In Figure 19, one of the 
interesting dynamic characteristics of the aircraft at this 
ang 1 e-of -attack (65 degrees) is the unstable, oscillatory 
Dutch-roll mode of the sailplane. This mode is signified by 
a pair of imaginary poles of the aircraft's characteristics 
equation oscillating at a frequency of 1.34 rad/sec. The 
near zero damping of Dutch-roll mode (.015), along with the 
low damping of the oscillatory roll-spiral mode (.91) might 
explain the observed flight oscillations. It is important 
to notice that the flight-determined dutch-roll oscillation 
frequency of 1.39 rad/sec is very close to the dutch-rol 1 
frequency of 1.34 rad/sec computed by DIGIKON, and plotted 
in Figure 19. One of the zeros of the of the 0/g a transfer 
function is shown to be in the right hand plane which is 
considered unconventional. But this characteristic does not 
appear to contribute to any. possible instability, because 


two S-plane root locus plots of 


and ff/ 
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as indicated by the gain versus real axis plot of Figure 
19, a very large pilot gain (about 3.5) is required to 
drive the roll-spiral mode unstable (pilot gain is'normally 
between .5 and 1). The Dutch-roll mode, on the other hand, 
has a tendency to become more stable as the pilot's gain is 
increase to a value of 1. A pilot gain of 1 can easily be 
achieved at such high angl e-of -at tack due to low aerodyna- 
mic damping and control surface effectiveness. 

Although the transfer function is considered to 

be less important than , a root locus analysis of it 

may help in presenting a clearer picture of the over-all 
dynamic response of the aircraft in the high angle-of- 
attack deep stall flight region. Figure 20 contains the S- 
p 1 ane root locus plot of /^ f for Flight 7. As indicated by 
Figure 20, the coupled rol 1-spiral mode has a tendency to 
become unstable if the pilot makes any effort to control 
the sideslip angle using his rudder. In this case any 
increase in the gain tends to stabilize the dutch-roll mode 
in a desirable manner. 
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CONCLUDING REMARKS 

A flight investigation of the controlled deep stall 
concept was conducted using a full scaled, modified, 
piloted SGS 1-36 sailplane, at ang 1 es-of -at tack of between 
-5 and 75 degrees. Despite the pioneering nature of this 
investigation, all of the primary objectives of the project 
were achieved, and for the first time a complete set of 
flight determined dynamic stability and control derivatives 
was obtained for an aircraft at such extreme angies-of- 
at tack . 

The aircraft was flown successfully in the deep 
stall region using large in-flight deflection of the all- 
movable horizontal stabilizer. This proved to be a viable 
technique for longitudinally controlling the aircraft in 
the deep stall flight regime. The NASA SGS 1-36 
demonstrated excel lent longitudinal and 1 atera 1 -d i rect i ona 1 
stability and control characteristics in the angle-of- 
attack range of 30 to 50 degrees, but displayed s ome 
undesirable 1 atera 1 -d i rect i ona 1 instabilities in the form 
of unstable Dutch-Roll oscillation above 55 degrees 
angle-of-attack. The magnitude of the oscillation appeared 
to be a function of angle-of-attack, and became more 
pronounced as the angle-of-attack increased. 
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In order to predict the dynamic behavior of the 
aircraft before the first deep stall flight, new techniques 
had to be developed for calculating the rotary derivatives 
at extreme ang 1 es -of -at tack. Post flight comparison of the 
predicted and the flight determined rotary derivatives, 
showed some of these techniques to be quite good while 
others proved inadequate. The comparison of static force 
and moment derivatives with that of the wind tunnel tests 
indicated some discrepancies in the low as well as the high 
ang 1 es -of -at tack regions, which could be attributed to the 
Reynolds number difference between the full scale and the 
quarter scale wind tunnel model. 

The research of the deep stall phenomenon is in its 
infancy. The present investigation has only scratched the 
surface in revealing the handling qualities and dynamic 
responses of aircraft in the deep stall flight regime. And 
more work is left to be done to gain more understanding of 
different phenomenon observed during the SGS 1-36 flights, 
the foil ow ing is a list of other possible areas were mo r e 
research can take place. 

1) A more detailed analysis of the data obtained from 
the present research could be conducted to provide more 
concrete evidence as to the causes of the unstable lateral - 
directional behavior exhibited by the aircraft at very high 
ang 1 es -of -attack. 

2) Detailed investigation of the separated flow 
through further test flights of the SGS 1-36 with the nose 
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boom removed to determine the effects of the vortesis shed 
by the nose bo om on the aircraft dyn am ics in deep stall 
regions. 

3) Spin testing of the SGS 1-36 in determining the 
possible spin recovery techniques using the in-flight 
deflection of the horizontal stabilator. 
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TABLE 1.- PHYSICAL CHARACTERISTICS OF THE 
NASA SGS 1-36 VEHICLE 


Total height, m (ft) 1.50 (4.92) 

Total length, m (ft) 7.27 (23.86) 

Wing - 

2 2 

Reference 6 actual planform area, m (ft ). 13.08 (140.80) 
Reference & actual span, m (ft) ...... 14.07 (46.17) 

Mean Aerodynamic Cord, m (ft) 0.97 (3.19) 

Reference cord, m (ft) 1.00 (3.28) 

Root cord, m (ft) 1.28 (4.20) 

Tip cord, m (ft) 0.58 (1.89) 

Taper ratio 0.451 

Aspect ratio 15.15 

Incidence at the root, deg 1 

Incidence at the tip, deg 0.065 

Sweep angle, deg 0 

Dihedral angle, deg 4 

Airfoil section at the root FX 61-163 

Airfoil section at the tip FX 60-126 

Horizontal tail - 

2 2 

Planform area, m (ft ) 1.23 (13.2) 

Span, m (ft) 2.41 (7.92) 

Average cord, m (ft) ^ . 0.51 (1.67) 

Root cord, m (ft) 0.56 (1.83) 

Aspect ratio • 4.76 

Leading edge sweep angle, deg 5 

Dihedral angle, deg ..... 0 

Stabilator range, deg 3 TED/76 TED 

Airfoil section NACA 64-012 

Vertical tail - 

2 2 

Planform area, m (ft ) 0.85 (9.10) 

Span, m (ft) 1.23 (4.04) 

Average cord, m (ft) 0.69 (2.25) 

Root cord, m (ft) 0.80 (2.64) 

Aspect ratio ........ 1.8 

Leading edge sweep angle, deg ....... 27 

Airfoil section NACA 64-012 
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TABLE 2.- PHYSICAL CHARACTERISTICS OP THE 
PRIMARY CONTROL SURFACES 


Ailerons - 

2 2 

Total planform area (both wings), m (ft ). 1.01 (10.90) 

Span, n (ft) 2.53 (8.30) 

Root cord, m (ft) 0.56 (1.84) 

Tip cord, m (ft) 0.39 (1.29) 

Average cord, m (ft) 0.40 (1.31) 

Hinge line, % average cord 72 

Range, deg . 12 Up/32 Down 

Elevator - 

2 2 

Total planform area, m (ft ) 0.42 (4.48) 

Root cord, m (ft) 0.23 (0.74) 

Tip cord, m (ft) 0.17 (0.56) 

Range, deg 20 Up/10 Down 

Rudder - 

2 2 

Planform area, m (ft ) 0.44 (4.72) 

Span, m (ft) 1.13 (3.70) 

Root cord, m (ft) 0.42 (1.39) 

Tip cord, m (ft) 0.30 (0.97) 

Average cord, m (ft) 0.39 (1.28) 

Hinge line sweep angle, deg 8 

Range, deg 30 Lft/30 Rgt 


TABLE 3.- MASS PROPERTIES OF THE NASA SGS 1-36 


Masses - 


Empty weight, N (lbs) 3002.4 (675.0) 

Useful load, N (lbs) 1089.8 (245.0) 

Maximum gross weight, N (lbs) 4092.2 (920.0) 

Inertias - 

Long, center of gravity, % MAC 33.4 

Measurement weight, N (lbs) 3883.1 (873.3) 

Ixx, Kg-m2 (slug-ft2) 1374.9 (1014.4) 

iYYr Kg-m2 (slug-ft2) 869.3 (641.3) 

Izz, Kg -m2 (slug-ft2) 2213.5 (1633.0) 

Ixz , Kg-m2 (slug-ft2) . 67.0 (49.4) 
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TABLE 4.- SELECTED INSTRUMENTATION RANGE AND ACCURACY 
FOR THE NASA SGS 1-36 VEHICLE 


Parameter Description 

Range 

Resolution 

(Bit/unit) 

Accuracy 

Angle of attack, deg 

-10 

9 

90 

0.1193 

0.5 

Angle of sideslip, deg 

-30 

9 

30 

0.0717 

0.2 

Pitch attitude, deg 

-42 

9 

42 

0.0822 

0.25 

Bank angle, deg 

-90 

9 

90 

0.1869 

0.54 

Elevator deflection, deg 

-25 

9 

25 

0.0439 

1.60 

Stabilizer deflection, deg 

0 

9 

100 

0.0597 

2.10 

Aileron deflection, deg 

-50 

9 

50 

0.0886 

4.60 

Rudder deflection, deg 

-50 

9 

50 

0.0606 

3.60 

Roll rate, deg/sec 

-60 

9 

60 

0.1171 

0.84 

Pitch rate, deg/sec 

-50 

9 

50 

0.0864 

1.30 

Yaw rate, deg/sec 

-50 

9 

50 

0.0870 

0.90 

Long, acceleration, g 

-0.5 

9 

0.5 

0.0012 

0.001 

Lateral acceleration, g 

-0.5 

9 

0.5 

0.0017 

0.002 

Normal acceleration, g 

-1.0 

9 

4.0 

0.0047 

0.002 

z 

Static pressure, hN/m 

587.9 

9 

936.7 

0.8700 

4.19 

Static pressure, lbs/ft 1232.0 

O 

9 

1963. 

0 1.8288 

8.77 

dm 

Dynamic pressure, hN/m 

0.0 

9 

16.7 

0.0200 

0.02 

Dynamic pressure, lbs/ft 

0.0 

9 

35.0 

0.0362 

0.04 

Static pres, rate, m/min - 

914.4 

9 

914.4 

1.8400 

1.83 

Static pres, rate, ft/min 

-3000 

9 

3000 

6.0460 

6.00 


TABLE 5.- MAXIMUM CRAMER- RAO BOUND TOLERANCES 
FOR PLOTTED DATA INCLUDING A SCALE 
FACTOR OF 3 


Parameter 

ITolerance 

CM„ 

0.002700 

CM q 

4.000000 


0.001200 


0.003000 

CN« 

0.013000 

cl o 

0.000100 


0.000400 

C1 «a 

0.000130 

Cn Ja 

0.000028 

C1 «r 

0.000022 

Cn Jr 

0.000018 

C1 P 

0.015000 

Cl, 

0.013000 

Crip 

0.015000 

Cn, 

0.010000 

a t 

0.000800 

CT «r 

0.000500 


0.000400 
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Figure 4: Stabilator control lever positioned for 60° stabiiator deflection. 
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Figure 5: SGS1-36 instrument panel with the automatic canopy ejection 
control handle. 





Figure 6: SGS1-36 nose cone with the boom installed. 
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Figure 5: SGS1-36 instrument panel with the automatic canopy ejection 
control handle. 
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Figure 7: The onboard data acquisition system. 



Figure 9: Basic controlled daap stall concept 
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Figure 11: NASA SGS1-36 deep stall flight program 





Figure 1 
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Figure 18: Continued. 
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Figure 18: Continued. 






Appendix A: This appendix presents the NASA SGS 1-36 non- 

dimensional aero-data package used for the real time 
simulation and initial linear analysis. The same 
data set was used as the initial starting values for 
the MMLE analysis. The static force and moment 
derivatives were obtained from the wind tunnel 
testing of the quarter scaled "free flight" model in 
stability axis and were transformed into body axis 
for presentation here. All the rotary aerodynamic 
stability derivatives are the predicted values as 
describe in the " Computed Aerodynamic Stability 
Derivatives" section of this report and are 


presented in body axis 
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CLB VS. ALPHA SGS 1-36 
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Appendix B: Linearized equations of motion as used by the Modi- 

fied Maximum Likelihood estimation program* 

The general fora of the equations is: 

*(«>*«> * A (t)x(t) + S(t)u(0 +S(f) 

r(t) * C (t)x(f) + D (f)u(f) + ff(t) *£(f)x(f) 
where A, 8, S, R r C, D f H# and E are defined by rel- 
ationships such as: 


A if (t) *AAf (/ (D XAHy+ALyW 

The a c and P c used are obtained from <* m and by the 
equations: 

« c « * aa^DCOFT) . £2$J£ 


where f 


KALP 

DCGFT 

XAN 

XAY 

ZAX 

ZAY 


Measured angle of attack 
Computed angle of attack 
Measured angle of sideslip 
Computed angle of sideslip 

Plow amplification factor for angle of attack 
Flight C*G. minus wind tunnel C*G. 

Longitudinal location of the normal acceleration 

sensor 

Longitudinal location of the lateral acceleration 

sensor 

Vertical location of the longitudinal acceleration 
sensor 

Vertical location of the lateral acceleration 


sensor 
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The nondimensional longitudinal matrices are: 
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The longitudinal dimensionalization matrices defined by the user routines are: 



CM - 


DM - 


1 

l 

"V 

1 


~ 1 
1 

l 

1 

i 

1 

n 

i ! 

1 

l 

1 


i 

1 

i 

i 

L j 

1 

1 

1 


si 

os 

qs 

as ! 

IL 

<f s c 
mgWS 

3* c 
"mg 2WJ 

1 


mg 

mg 

mg 

mg 

mg 

.51 

mg 

1 

1 

-51 
_ mg 

•51 

mg 

as 

"mg 

-SI 

mgj 


EM - 


H.U - 





fl 1 


1 

1 

1 

17 


i i i 

1 * 

i i 
* ! 

1 

1 

1 

i 


i 1 i 

1 

i 

1 

1 

1 

1 


j 1 & 


j 

as 

cs 

cs 

qs ! 


1 

J 

mg 

mg 

mg 

mg i 


l & 

‘ 1 

i 

! 

os 

- 

as 

'mg 

as 

'mg 

3 0 

<^[co 

1 
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The FAI matrix is filled with l*s . 


AL - 

1 ^(-cos co sin 0 cos a 
- cos 0 sin <x^ 

0 0 

cos <p 0 

The BL , SL, CL , DL , HL . and EL matrices are filled with 0's . 


0 

0 

0 


are: 


Lateral-directional . - The nondimensional lateral-directional matrices 


AN - 


BN - 








r 


0 


SN - 



C 


n 


0 


L 
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The dimensionalization matrices defined by the standard aircraft routines for a 


lateral-directional case are: 


AM - 


■ft* 

b 

mV W 

7 s b 
mV IV 

l x 

Tfsb b 

17 IV 

?sb b 
I IV 

X 


5 ab b 

rw 

Tsb b 

T7 W 

. 1 

1 

1 


AM - 




" in. 
mr 

CS -p 


JJSf' 

l x 

*X 


J x 

* 


csb 

v* 


1 

1 

1 

1 





CAf - DM - 


“ 

1 

1 

* 


’ 1 

1 

1 

1 * 

1 

'7 

V 

1 


e 

1 

1 

i 

1 

1 

t 

A 

1 

1 


1 

a 

1 

1 

1 

1 

l 

1 


1 

1 

1 

1 

1 

1 

1 

1 





r 



cs 

2l 

OS 

££ 

Sl 

7$ b 

7s b 



mg 

mg 

mg 

mg i 

mg 

mg IV3F 

mg KS 

1 1 


1 

1 

i 

1 

i j 

j 1 

1 

t 

1 


i_ 

1 

1 

, ; 

1 





X 


*■ - 

i 1 

1 

1 

lj 
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HM - EM -r 



The FM matrix is filled with l's. 

AL - 

0 sin a -cos a cos cp cos 9 
0 0 0 0 

0 0 0 0 

0 1 cos <p ten 9 0 

The BL , SL , RL , CL . DL , HL , £L , and FL matrices are filled with 0’s . 



Appendix C: The following equation was developed to estimate 


the value of Cl for the angles of attack of 
P 

between 30 and 90 degrees. 

It is assumed that since the wings are complete* 
ly stalled above 20 degrees angle of attack, it 
is reasonable to consider the normal force (or 
drag at high angles) to be the primary damping 
force. 


- b /2 




dy 



S 



We have : 


C * C - 2y/b(C -C ) for +y (1) 
local r r t 


Let: 


X-2/b(C -C ) (2) 

r t 


And rolling moment Cl to be: 


r /2 

Cl - JCN . y . q • C dy/q.s.b (3) 

-b^ local local ®* 
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For positive p, 

v » P y + V_ 

local 

q « 1/2 p (P y + v_) 2 

local 


(4) 

(5) 


and we know that: 


2 

q„ - 1/2 p 


cn - f ( a ) 
local 

a » a„ + 

local 

da- f ( p y ) 


( 6 ) 

(7) 

(8) 


substituting the equations 1, 5, and 6 into 
equation 3 we have: 


Cl - CN 


{ y"yd > y+v.) 2 (c -\ y ) dy + 

0 0 r 2 

fy{Py+v m uc + X y ) dy j/ v„ s b (10) 


integrating equation 10 we will have, 

2 3 

Cl « CN/ v„ s ( b /6 C P - b /16 X P ) (11) 

r 

and finally. 


2 3 

Cl ■ CN/2 s b ( b /$ C - b /16 X ) (12) 

P r 
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Sine* the value of CN as a function of a has 
been determined from the wind tunnel data the 


parameter Cl can now be calculated using 
P 

equation 12. 


